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To create a stable signal from a cryogenic sapphire maser frequency standard, the 
frequency-temperature dependence of the supporting Whispering Gallery mode must 
be annulled. We report the ability to control this dependence by manipulating the 
paramagnetic susceptibility of Fe 3+ ions in the sapphire lattice. We show that the 
maser signal depends on other Whispering Gallery modes tuned to the pump signal 
near 31 GHz, and the annulment point can be controlled to exist between 5 to 10 
K depending on the Fe 3+ ion concentration and the frequency of the pump. This 
level of control has not been achieved previously, and will allow improvements in the 
stability of such devices. 
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I. INTRODUCTION 



The Whispering Gallery Mode Sapphire Maser is a stable microwave oscillator, which 
makes use of paramagnetic Fe 3+ ions in ultra-low loss HEMEX cylindrical sapphire crystal 
at low temperature, and has been described in detail previously^-. The fractional fre- 
quency instability has been demonstrated to be as low as cr y (ls < r < 100s) = 10~ 14 . This 
instability is normally measured at the frequency-temperature turnover point (or annulment 
temperature) of a high Q-factor (> 10 9 ) Whispering Gallery Mode (WGM) excited in the 
sapphire, where the effects of temperature fluctuations on frequency are nullified to first 
order by the magnetic susceptibility of residual paramagnetic impurities^. Many publica- 
tions describe the application of this self-compensation technique^— , and all show that the 
annulment temperature is dependent on the presence and relative concentrations of param- 
agnetic ions such as Cr 3+ with an Electron Spin Resonance (ESR) at 11.4 GHz, Mo 3+ (100 
GHz) and Ti 3+ (1 THz). The ions are substitutionally included in the crystal lattice during 
the crystal growing process and occur unintentionally. 

In this work we show that the frequency-temperature annulment for the Maser depends 
predominately on the Fe 3+ ions and can be controlled by manipulating the number of ions 
involved in the maser process. We compare the behavior of two crystals - CI (concentration 
of active Fe 3+ ions = 10 ppb) and C2 (100 ppb). We also report the observation of an upper 
limit in temperature of 30 K for operation of the maser. This limit is explained using a 
Boltzmann distribution for the active ion, and we show that the populations of the energy 
levels are so close at this temperature that the effect of the pump at 31 GHz used to excite 
the maser at 12 GHz is negligible. 

Controlling spins in such high-Q cavities could also have other potential applications. For 
example, it has recently been suggested that HEMEX sapphire with paramagnetic impurities 
could have applications for quantum measurement at millikelvin temperatures^ 1 ^, such as 
qubits or quantum memory applications^— . 

II. MASER DESCRIPTION 

The maser scheme is based on the three energy levels of Fe 3+ ions at zero applied DC 
magnetic field in the sapphire lattice of a WGM resonator (Fig. [Q). The concentration of 
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active ions used to create the effect is about 10-100 ppb, whereas the total concentration 
is on the order of 1 ppm. The bandwidth of the ESR is A^ Fc 3+ « 27 MHz, where the WG 
mode width is only Ai/ WG ~ 10 Hz at 4.2 K. 




FIG. 1. Principle of operation: We apply a signal at 31.34 GHz, coinciding with a WG mode, to 
pump the ions from the lowest energy level |l/2) to the third 1 5/2} . The ions then (non-radiatively) 
relax to the second energy level |3/2). A population inversion is obtained between the two lowest 
energy levels, and stimulated emission can be achieved at 12.04 GHz, which is enhanced by another 
WG mode coincident in frequency. 



The maser signal frequency u op is fixed at the frequency z/^g of the WG mode involved 
in the process. It is well known that the frequency %q is strongly sensitive to temperature 
variations and needs to be precisely controlled. In this kind of resonator, iavg possesses a 
frequency-temperature turnover point (also referred as the annulment or annulment temper- 
ature) above 4.2 K, although annulment temperatures as low as 90 mK have been recently 
measured^. The frequency variations can be described by^: 

UWG ~ Uo = AT* + C(T, u) (1) 

where vq is the mode frequency at K without ions, and AT 4 is the thermal sensitivity due 
to the sapphire dilation (A « 10~ 12 K -4 for our modes). The other term C(T, v) is thermal 
sensitivity due to the paramagnetic ions and is defined by: 

C(T, v) (2.rf)V^-,) 

This term depends in reality on both temperature and operational frequency. The annulment 
temperature is then calculated by considering the maxima of Eqn. [21 where the index i labels 
the different paramagnetic species present in the lattice. The magnetic filling factor rji of the 



mode describes its magnetic field distribution in relation to the i th species of paramagnetic 
ion, and Xo is the DC magnetic susceptibility of that ion. As the operating frequency is 
fixed at around 12.04 GHz, only the effect of the Fe 3+ is considered (i.e. the summation in 
Eq. [2] is only performed over a single ion i = Fe 3+ ), and the annulment temperature (T inv ) 
can then be expressed as follows^ 3 -: 
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where the Bohr magneton \xb = 9.274 x 10~ 24 J.T -1 , the Boltzmann constant ks = 
1.38 x 10 -23 J.K -1 , the Lande g-factor g « 2, a\ 2 ~ 2 is a characteristic constant of the 
absorption rate of the ions at 12.04 GHz (calculated from the Fe 3+ Sapphire spin Hamiltonian 
parameters)^ 3 - - — , the spin-spin relaxation time r 2 is related to the width at half resonance 
of the ESR (Au 12 = 1 /ttt 2 ). If two spins are in phase, and one absorbs a photon causing a 
phase shift, t 2 is the necessary time for the spins to be in phase again. 
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FIG. 2. The active maser bandwidth over which we can generate a signal is of order 180-190 MHz 
for crystal C2. Masing occurs at discrete frequencies where WGMs exist between 11.95 to 12.15 
GHz, which are in turn pumped at discrete frequencies between 31.24 to 31.44 GHz where higher 
frequency WGM's exist. 



Fig. |2] shows the bandwidth for active maser operation in crystal C2, which is more than 
186 MHz for the 12 GHz and 31 GHz ESR levels. We can note from Fig. [2] that the maser 



can operate in a multimode configuration. For example, when the pump frequency is in the 
range 31.350 GHz - 31.375 GHz, two signals can operate at once at 12.037 GHz and 12.149 
GHz. In the rest of the paper, the characterization of the turnover temperature has been 
investigated for single-mode maser operation. We consider then in this paper Az/ 12 = 27 
MHz, the width measured by Bogle and Symmons 2 ^. N int is the concentration of the active 
ions involved in the annulment temperature, in units of ion/m 3 . However, it is convenient 
to express it as a fraction of the AI2O3 in the crystal, i. e. 2.35 x 10 22 ion/m 3 = 1 ppm. 
This concentration is related to the total ion concentration in the lattice as follow: 

N int = ( x N total (4) 

where ( is the fraction of the ion concentration involved in the proces (0 < ( < 1) and N to tai 
the total ion concentration in the lattice. When £ = 0, no frequency-temperature annulment 
is observed, and when £ = 1 all the ions in the lattice participate in the annulment process. 
Fig. [3] shows that the annulment temperature T; ni) increases with the total number of ions 
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FIG. 3. Evolution of the annulment temperature with the total Fe 3+ concentration for different 
values of the fraction 

N total ■ 

III. MEASUREMENTS 

To characterize the maser signal, we amplify it by 30 dB and then compare it to a signal 
from a synthesizer referenced to a hydrogen maser. The beat note is then amplified and 
band pass filtered before being counted (Fig. [4]). The accuracy of the temperature control 
at the resonator cavity is about 2 mK. 
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FIG. 4. Experimental setup: The crystal is a cylindrical HEMEX sapphire, 30mm height and 50mm 
diameter, placed in a cylindrical cavity and cooled to 4-2 K using a cryocooler. 

A. Behavior of CI with low concentration of Fe 3+ 

The maser characterized in this section oscillates at the WGH17 o,o mode frequency (12.038 
GHz) of crystal CI, and can be excited using a pump signal corresponding to any of several 
WG modes around 31.34 GHz. The frequency of the WG mode that enhances the maser 
signal shows a turnover temperature T inv around 7.8 K, with no dependence on the input 
power applied to the crystal. The evolution of the maser frequency with temperature is 
shown in Fig. [5j 
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FIG. 5. Evolution of the maser frequency (for CI) at 12.038 GHz with temperature for a fixed 
pump at 31.339,734 GHz. From the fit, the annulment temperature is 8.193 K. 

The excitation frequency of each pump WG mode was varied, and the evolution of the 
maser frequency with temperature was observed at various levels of pump power. The 
frequency was changed in steps of 250 Hz at various levels of pump power. As Fig. [6] shows, 
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Ti nv is independent of the pump power and the pump frequency for each pump mode. For 
example, the 31.339 GHz pump mode has an annulment temperature which stays constant 
at 8.25 K. 
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FIG. 6. Evolution of the annulment temperature for different values of the pump power and fre- 
quency. The pump mode chosen was at 31.339 GHz. 

Despite the fact that T inv is independent of the input power and frequency of the pump 
signal (for a particular choice of pump mode) , its value does change when a different pump 
WG mode is selected, as shown in Fig. [71 The Ti nv point is measured for the pump frequency 
range 31.320 - 31.339 GHz. No turnover point is observed at higher frequencies, due to the 
annulment temperature residing outside the measurable temperature range (below 4 K). T inv 
can be calculated for each pump by a simple second order polynomial equation, as shown 
in Fig. [5j 

For each pump at 31 GHz, the population of active ions N maser used to create the maser 
signal is different from that fraction of ions that determine the Ti nv for the mode. The 
active ion concentration was deduced for each T inv from the maser output power calculated 
by Benmessai et al.— , as shown by Eq. [51 
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The concentrations N maser and N int are the same only if the ESR is saturated. N n 
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FIG. 7. Evolution of the maser signal frequency with temperature for different pumps at 31.3xx 
GHz. For each pump T{ nv is independent from its frequency and the power. 

is calculated from the rate equations and the Boltzmann law at saturation, and N int is 
calculated from WG mode behaviour. The ion inversion ratio I ra uo is defined as follows: 

(AjV 23 - Aiy 12 )(Aiy 13 + AjV 23 ) 

AiV 12 (5AiV 1 3 + 2AiV 2 3) [) 

where the different AiVy are the normalized population difference between the energy levels 
% and j and are calculated as shown by Benmessai et al.— . The relaxation time T\ ~ 10 ms 
at 8 K is the spin- lattice relaxation time. This time is considered as constant around T inv . 
In reality it undergoes a direct relaxation process and follows T" 1 when T < 15 K— ~— . 
Here, (3 out is the coupling of the mode at 12.04 GHz, Q is its unloaded Q-factor, V e ff is the 
effective mode volume, t] is the filling factor. The mode parameters are constant over the 
temperature range in which we characterize the signal, and are summarised in Table [B 



TABLE I. Measured Parameters for the WGH170 mode at 8 K. The maser signal was operating 
on the lower doublet of the mode. 



Mode 




v WG (GHz) 


Pout 


Qo x 10 6 


Veff (m 3 ) 


V 


WGHi^co 


7.8 


12.038136 
12.038137 


0.200 
0.009 


600 


lQ- r ° 


0.99268 
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By definition, £ is the ratio between the power in the unsaturated regime and saturation 
(0 < £ < 1). Note that for the pump at 31.337 GHz where the maser is at saturation, £ 
should be equal to unity. From the results shown in Fig. [71 where T inv is different from one 
pump to another, it is possible to calculate the corresponding concentration for each T inv 
from Eq. (121). It is also possible to calculate the same concentration from Eq. (j5J). Knowing 
the maser output power, the result are shown in Table [Til 



TABLE II. Summary of the calculated concentrations 



Vpurnp (GHz) 


Ti n u (K) 


Power (dBm)^ 


N int (ppb)i 


c 


Nmaser (ppb)£ 


e 


31.320 


9.834 


-62.59 


13.46 


0.0135 


7.74 


0.6180 


31.337 


8.656 


-60.50 


7.11 


0.0071 


11.09 


1.0000 


31.339 


8.169 


-65.50 


5.32 


0.0053 


3.57 


0.3162 


31.347 




-69.00 d 






1.90 


0.1412 



a Measured at the output of resonator 
b Calculated from Eq. [3] 
c Calculated from Eq. [5] 
d at 9 K 

Eq. (jSJ) is evaluated when the maser signal is maximum (i.e. at saturation). As the maser 
signal is oscillating at one mode (WGH^o.o) for all the pumps, it is clear that for each pump 
the amount of active ions is different, leading to a calculation of different concentrations. 
Each pump has a different configuration of coupling and field distribution in the crystal. So, 
each pump will interact with a different quantity of ions, and whether saturation is reached 
is determined by the coupling. 

B. Behavior of C2 with high concentration of Fe 3+ 

It has been shown previously^ that the concentration of Fe 3+ ions in sapphire can be 
increased significantly by annealing in air, causing conversion of Fe 2+ impurities to Fe 3+ . 
Originally the two crystals CI and C2 exhibited similar properties of very low concentration 
of Fe 3+ impurities. However, C2 went through a series of additional annealing and was 
thus transformed from the low concentration regime (as reported here for CI) to the high 
concentration regime. The improvement in ion concentration is permanent, and once in the 
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high concentration regime the resonator will exhibit the effects reported in this paper. Thus, 
annealing a resonator is beneficial in two respects; due to the increased active ion population 
it will exhibit a higher maser output power, leading to a reduction in the Schawlow-Townes 
thermal noise limits, as well as leading to a potential improvement in the curvature of the 
frequency-temperature dependence as is reported in this section. 

For the crystal C2, the maser operation is very different as the concentration of the active 
ions is higher, resulting in a higher maser output power of —40 dBm, compared to —56 dBm 
for CI. In addition, maser signals are observed not only at the WGHi 7i o,o frequency, but 
also at different modes whose frequencies are within the ESR bandwidth^. Like CI, all the 
pump modes for which a maser signal can be excited in C2 have a temperature turnover 
point independent of the applied power. Characterizing the turnover point for this crystal 
reveals more complex effects than CI. 
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FIG. 8. Evolution of the maser signal frequency at 12.037 GHz with temperature for different value 
of the same pump at offsets (in Hz) from 31.349,780,000 GHz. 

Figure [8] shows the signal at 12.037 GHz oscillating for a pump applied at 31.349 GHz. 
Each curve represents the maser frequency dependence as a function of temperature for 
different pump frequencies. Here, the maser frequency shows a strong dependence on the 
pump tuning. We note that when the frequency u pump < 31.349,782,315 GHz, there are 
three Ti nv - one near 6 K, the second around 8 K, and the third at 10.5 K. When v pump > 
31.349,782,415 GHz, there is only one T inv around 8 K and the system behaves in the same 
classical fashion as the first crystal. 
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We can also divide the figure into four zones, with each zone defining a turnover temper- 
ature. From the results shown in Fig. HI we can draw the evolution of the different turning 
point with the pump frequency (cf. Fig. [H]). 
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FIG. 9. Evolution of the turnover temperature at 12.031 GHz with different values of the same 
pump at offsets (in Hz) from 31.349,780,000 GHz. 

The annulment exists for multiple temperatures only for the crystal C2 with high Fe 3+ 
concentration. This effect appears to be nonlinear, showing that a variation in the active 
ion population is induced as the temperature is varied, which inturn changes the turning 
point temperature as the temperature is varied. 

The annulment temperature depends on the distribution of the ions and the population 
difference between their levels. When the annulment temperature for the WGH 17i0i o mode is 
measured using a network analyzer, no change is seen due to the small population difference 
at this frequency (5% of the total active population). However, when the maser is operating, 
i.e. being pumped, the configuration of the ions is completely different: population inversion 
exists between the 1 1/2) and |3/2) levels, and saturation is seen at 31 GHz where the pop- 
ulation of the states |l/2) and 1 5/2) are equal. The population inversion rate is determined 
by the choice of pump mode near 31 GHz. For each pump mode, the strength of the maser 
signal depends on its coupling and spatial field distribution in the crystal. Strong coupling 
allows an optimum exchange of energy with the ions, and field distributions occupying more 
volume in the crystal increase the interactions with different classes of ions. Considering 
the inhomogeneous broadening of the line at 31 GHz in addition to these factors, we can 
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understand why Ti nv is different from one pump to another. The Fe 3+ spins act as individual 
packets selected by the pump WG mode corresponding. Thus, the pump signal selects a 
different quantity of ions among the broadened ESR frequency. 

For C2, the case is more complex because T; n „ is not only different from one pump 
to another for the same maser signal, but is different for the same pump mode when the 
frequency is tuned. The maser frequency shows more than one temperature turnover point 
for some pump frequencies, and no annulment for others. In the case where there is more 
than one, the curvature of annulment is significantly reduced, which has the potential to 
improve the stability and reduce the requirements of temperature control of the maser. 

Optimum Operation 

The characterization of the maser frequency with temperature was performed at a fixed 
pump frequency. However, the pump WG mode is temperature dependent as well as the WG 
modes near 12.04 GHz. To optimize the maser operation, it is necessary to tune the pump 
frequency for each operating temperature. This optimum frequency corresponds to maxima 
or minima of the maser frequency depending on the Tj m , being measured. Fig. HO] shows 
how the maser frequency varies with temperature for different pump modes. The signal at 

120600 , , , , , , , , , 




4 5 6 7 8 9 10 11 12 
Temperature (K) 

FIG. 10. Evolution of the maser signal frequency with temperature for different pumps. 
12.014 GHz has a turnover point at 10 K, the one at 12.029 GHz is at 10.3 K, and the one 
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at 12.037 GHz shows no turnover temperature at all. The measurements were performed at 
different power levels, which were found to have no effect on the turnover temperature. 



IV. MASER TEMPERATURE OPERATION LIMIT 

In order to characterize the upper temperature limit T opera ti on for maser operation, we 
adjust the pump frequency for each measurement as per the last section. The results are 
summarised in Table HTT1 The 12.037 GHz maser signal oscillates until 28 K when pumped 
at 31.349 GHz, and 30 K when pumped at 31.352 GHz. The power of the maser signal is 
observed to decrease as the temperature is raised. This is due to the fact that the populations 
of the Fe 3+ energy levels follow a Boltzmann distribution at thermal equilibrium. At higher 
temperatures, the difference in population between levels is reduced, resulting in less efficient 
masing. At high temperature (20-30 K) the number of ions pumped at 31.3 GHz is lower 
than at 4.2 K, leading to a smaller population inversion ratio. This directly affects the maser 
power and causes it to drop as temperature rises. In addition, the absorption at 31.3 GHz 
is a result of a mixing between the 1 1/2) and 1 5/2) spin states of the Fe 3+ ion, which occurs 
with very small probability^' 3 ^' 3 ^. The population difference at 31.3 GHz is high enough 
to excite the maser signal with this probability at low temperatures (<30 K), but not at 
high temperatures (>30 K). At higher temperatures still, characterization of the WG modes 

TABLE III. Different annulment and operation limit temperatures for optimal operation of the 
maser 

Vpump (GHz) V maser (GHz) Tinv (-K) ^operation (K) 

31.310 12.014 9.407 15 

31.319 12.029 9.910 26 

31.334 12.029 8.140 23 

31.349 12.037 - 28 

31.352 12.037 - 30 



around 12.04 GHz showed absorption effects up to 140 K. The strength of the absorption 
decreases when the temperature is increased, however no maser operation has been observed 
at these temperatures. 
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CONCLUSION 



We have demonstrated how Fe 3+ ions in a sapphire maser influence the frequency- 
temperature turnover point of the maser signal. We compared the behaviour of two crystals, 
one with a low active ion concentration and the other with a high concentration. The first 
crystal showed classical behaviour where Ti nv is independent of the pump power and fre- 
quency, but was different between pump modes due to different packets of ions being excited 
for each pump. The population of ions that creates the maser causes a change in the intrin- 
sic annulment temperature Ti nv for the mode. The crystal with higher Fe 3+ concentration 
showed a more complex behaviour, where T inv was different not only from one pump mode 
to another, but showed many turnover points for some choices of pump, and in some cases 
none at all. Thus, the curvature of the annulment point in some cases could be reduced, 
which would also reduce the requirements for temperature control when generating a stable 
frequency. 
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